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(3J[) Abstract 



Let $ : 1 -> R be an arbitrary continuously differentiable deter- 
ministic function such that |$| + |$'| is bounded by a polynomial. 
In this article we consider the class of stochastic volatility models in 
which {Z(i)}te[o,i]; the logarithm of the price process, is of the form 
Z(t) = L $(X(s)) dW(s), where {^(s)} s e[o,i] denotes an arbitrary 
centered Gaussian process whose trajectories are, with probability 1, 
Holder continuous functions of an arbitrary order a € (1/2,1], and 
where {T / ^(s)} sS [o.i] is a standard Brownian motion independent on 
{^(s)} s£ [0,i]- 

First we show that the critical Holder regularity of a typical tra- 
jectory of {Z(t)} te [ 01 ] is equal to 1/2. Next we provide for such a 
trajectory an expression as a random series which converges at a ge- 
ometric rate in any Holder space of an arbitrary order 7 < 1/2; this 
expression is obtained through the expansion of the random function 
s n- &(X(s)) on the Haar basis. Finally, thanks to it, we give an 
efficient iterative simulation method for {Z(t)} t£ ^ 01 ]. 
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1 Introduction 

Stochastic volatility models are extensions of the well-known Black and Sc- 
holes model. Hull and White [15] and other authors working in the field 
of mathematical finance (see for instance [22] and [E]), introduced them 
starting in the mid-80s, in order to account for volatility effects of exoge- 
nous arrivals of information on the price of an underlying asset. Later, for 
making such models to be more realistic, it has been proposed (see for ex- 
ample \TT\ [131 [TH 121] ) to replace the Brownian motion governing the 
volatility in them, by a more flexible stochastic process. That is why, in this 
article we consider the general class of stochastic volatility models, in which 
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the stochastic process {^(i)}te[o,i] denoting the logarithm of the price of an 
underlying asset, is defined for each time t G [0, 1], as, 

Z(t)= f\(X(s))dW(s), (1.1) 
J o 

where 

• {^(s)} s£ [o,i] denotes an arbitrary real-valued centered Gaussian pro- 
cess whose trajectories are, with probability 1, Holder continuous func- 
tions of an arbitrary order a G (1/2,1] i.e. there is a positive and 
finite random variable C such that one has for almost all u and for all 
si,s 2 G [0, 1], 

\X( 8l ,u)-X(82,u)\ < C(a;)|si-S2| a ; (1.2) 

it is worth noticing that, thanks to the Gaussianity of {-X"(s)} s e[o,i] , the 
random variable C can be chosen in such a way that all its moments 
are finite. 

• $ is an arbitrary deterministic real-valued continuously differentiable 
function on the real line, which satisfies as well as its derivative <&', for 
all x£R, 

\<&{x)\ + |$'(z)| < c(l + |x|) L , (1.3) 
where c > and L > are two constants (non depending) on x. 

• {W(s)} S £[o,ii denotes a standard Brownian motion independent on the 
process {-X"(s)} S £[o,il- Observe that thanks to these assumptions, the 
stochastic integral in (jl.ip is well-defined, moreover, conditional on 
the cr-algebra, 

g x = v{x(8),o<8<i), (i.4) 

{•Z(i)}te[o,il is a centered Gaussian process with a covariance function 
given, for all t\,t2 G [0, 1], by, 

/•min(ti,t 2 ) 

E(Z(t 1 )Z(t 2 )\Gx)= \$(X( S ))\ 2 ds. (1.5) 

■/ 

Observe that the models in (jl.ip are generalizations of the multifractional 
stochastic volatility models introduced in [2j , which in turn extend the frac- 
tional stochastic volatility models studied in [13\ I14| . Recall that a mul- 
tifractional stochastic volatility model is defined through (jl.ip in which 
{X(s)} sg [o j i] denotes a multifractional Brownian motion [6l [20], in other 
words, roughly speaking, a fractional Brownian motion with a smooth time- 
varying Hurst parameter H(s); a survey on the latter Gaussian process can 
be found in [I] for instance, also we refer to [TJ El [9] for some interesting 
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applications of it in the field of finance. Recall that the main feature of 
a multifractional stochastic volatility model, is that the local roughness of 
the volatility process {5>(JT(s)} s6 r in other words its pointwise Holder 
exponent (see (|2.ip for a definition of this exponent), can be prescribed via 
the Hurst functional parameter H (•) and thus is allowed to change over time 
(that is from point to point). 

The remaining of the paper is structured in the following way. In sec- 
tion[2l we show that there exists a modification of {#(t)}te[o 11, also denoted 
by [o,i] j whose trajectories belong, with probability 1, to any Holder 

space C 7 ([0, 1]) of an arbitrary order 7 < 1/2; on the other hand, under the 
assumption that <I> vanishes only on a Lebesgue negligible set, we show that 
the pointwise Holder exponent of {-Z(i)}te[o,i]> at any point to G [0, 1] such 
that E|X(to)| 2 > 0, is almost surely equal to 1/2; observe that, under the 
same assumption on $, the latter result implies that, with probability 1, 
the trajectories of {Z(t)} te [ ,i] fail to belong to C 7 ([0, 1]) when 7 > 1/2 and 
{X(s)} sg [ 0i i] is not an almost surely vanishing process (i.e. we do not have 
for all s G [0,1], almost surely, X(s) = 0). In section [31 by expanding the 
random function s 1— > &(X(s)) on the Haar basis of L 2 ([0, 1]), we introduce 
a random series representation of {Z(t)} tg j ,i] j f° r which the convergence 
holds at a geometric rate, almost surely in C 7 ([0, 1]), for all 7 < 1/2. Fi- 
nally, thanks to the latter nice representation of {Z(t)} t6 j ,i] 1 we gi ye i n 
section 0] an iterative algorithm which allows to efficiently simulate this pro- 
cess. 



2 Holder regularity of the log price process {Z(t)} te [Q^ 

Let us first recall the definition of a Holder space of order 7 G [0, 1]. 

Definition 2.1 For any 7 G [0, 1], the Holder space C 7 ([0, 1]) is defined as 
the Banach space of the continuous real-valued functions u which satisfy, 



Mil) - u(t 2 )\ 
sup < +00. 



It is equipped with the norm, 



\u(ti) — ^(^2)! 

Mlc7([o,i]) = IMloo + sup — — — , (2.1) 



o<ti<« 2 <i ri — *2| 7 



where \\u\\oo = sup fe [ 01 ] \u(t)\. 



Observe that C°([0, 1]) reduces to C([0, 1]) the usual Banach space of 
real- valued continuous functions over [0, 1]. The main goal of this section is 
to prove the following two theorems. 
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Theorem 2.1 Let {^(i)}te[o,i] be ^ e stochastic process defined in \l. 
then there exists a modification o/{Z(i)} fg [ 0i i], a/so denoted by {Z(t)} te \o,i]> 
such that, with probability 1, for all 7 < 1/2, i/te trajectories of {Z(t)} tE \o,l] 
belong to the Holder space C 7 ([0, 1]). 

Recall that when a trajectory i 1— )■ Y" (i, w) of an arbitrary real- valued stochas- 
tic process {Y(t)}te[o,i] is a continuous and non-differentiable function at a 
point to £ [0, 1], then py (io> w), the pointwise Holder exponent of 1 1— >■ Y~(t, cj) 
at to, is defined as, 

x / c rn 11 r r(to + /t,a;)-y(to,a;) 1 

py(to,w) = sup< p € [0,1] : limsup — = f- 

[ t +Ae[o,i], h-+o \ n v J 

(2.2) 

Theorem 2.2 Let {^(£)}te[0,i] be the modification of the process {Z{t)} t& \p,i} 
introduced in Theorem \2.1\ and let pz be the corresponding pointwise Holder 
exponent; then, under the assumption that $ vanishes only on a Lebesgue 
negligible set, one has, 

¥(p z (t ) = 1/2) = 1, V t € [0, 1] s.t. E \X(t )\ 2 > 0. (2.3) 

Observe that a straightforward consequence of Theorem 12. II is that, one 
has for all to 6 [0, 1], 

P(pjj(to) ^ V2) = 1- (2-4) 

Therefore, for obtaining Theorem 12.2^ it is sufficient to show that Theo- 
rem 12.11 holds and that, under the assumption that vanishes only on a 
Lebesgue negligible set, one has, 

P(pz(«o) < 1/2) = 1, V t € [0, 1] s.t. E \X(t )\ 2 > 0. (2.5) 

The proof of Theorem 12.11 mainly relies on the following two lemmas. Ob- 
serve that the first one whose proof can be found in |16| for instance, is a 
refined version of the usual Kolmogorov criterion allowing to show that a 
stochastic process has a modification with almost surely continuous trajec- 
tories; also observe that the second one, which can easily be proved, is a 
standard result on centered real-valued Gaussian random variables. 

Lemma 2.3 (Kolmogorov-Centsov criterion) Assume that the positive 
real number T is arbitrary and fixed. Let {M (i)}te[o,Tl be an arbitrary 
stochastic process such that for all (si,«2) £ [0,T] 2 , one has, 

E\M( Sl ) - M(s 2 )| T < c|si - s 2 | 1+/3 , 

where c, r and f3 are three positive constants. Then there exists a modifi- 
cation of {M(t)} te <o,T] denoted by {M(£)}te[o,T]j whose trajectories are with 
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probability 1, Holder functions of any arbitrary order 7 £ [0, /3/t); in other 
words, 

/ gup \mm-m( S2 )\ < \ =1 

yO<si<s 2 <T 1^1 — S2I 7 J 

Lemma 2.4 (equivalence of Gaussian moments) Let G be an arbitrary 
centered real-valued Gaussian random variable. Then for all positive real 
number t, one has, 

E|Gf = c T (E|G| 2 ) T/2 , (2.6) 
where c T denotes the positive constant, only depending on t, defined as, 

2 r/2 T (r±l) 



r being the usual "Gamma function", defined as T(n) = Jq + °° x u 1 e x dx, 
for all u £ (0, +00). 

Proof of Theorem I2.lt Let 7 G (0, 1/2) be arbitrary and fixed, we 
assume that the positive real number r has been chosen in such a way that 

7 < \ - i- (2.7) 

Recall that the cr-algebra Qx has been defined in (|1.4p and that conditionally 
on it, the stochastic process {^(^)}te[o,i] has a centered Gaussian distribu- 
tion with a covariance function given by (|1.5|) . Therefore, one has, for any 
real number r > and for each t\, i 2 G [0, 1], 

E|Z(ti)-Z(t 2 )| T = E{E(|Z(*i)-Z(t 2 )| T |fe)} 

= dEKEdZ^)-^)! 2 !^))" 72 }, (2.8) 

where c\ > is the constant c T in (|2.6p . moreover, 



Edz^o-z^)! 2 !^) 

= E (IZ^)! 2 !^) + E (|Z(t 2 )| 2 |fo) - 2E (z(h)z(t 2 )\g x ) 

= / |$(X(s))| 2 ds+ / \<S>(X(s))\ 2 ds-2 / |$(X(s))| 2 ds 
io io 

/•max(ii,t 2 ) 

= / |d>(^( S ))| 2 d S . (2.9) 

J min(ti,t2) 

Next, in order to bound E (\Z(ti) — Z(t2)\ 2 \Gx) , one sets 

11X1100= sup \X(s)\; (2.10) 
te[o,i] 
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observe that ||^||oo is an almost surely finite random variable, since the 
trajectories of {^(s)} s e[o,i] are > with probability 1 continuous functions. 
The Gaussianity of the latter process as well as the almost sure finiteness of 
|| X | |oo ; imply that (see [T7j for instance), all the moments of this random 
variable are finite as well, namely, for each real number p > 0, one has, 

E(||X[p <oo. (2.11) 

Next, one sets, 

Ci= sup |*(af)|, (2-12) 

observe that C\ is an almost surely finite random variable since $ is a 
continuous function, moreover, in view of (|1.3p and (|2,lip all its moments 
are finite, namely, for each real number q > 0, one has, 

E (Cf) < oo. (2.13) 

Next, combining (j2?9j) with (j27i"0j) and (j27f2j) . one gets that, 

E (\Z(h) - Z{t 2 )\ 2 \g x ) < Cf |ti - t 2 \, (2.14) 

therefore (|2.8p entails that 

E \Z(h) - Z(t 2 )\ T < c 2 \h - t 2 \ T ' 2 , (2.15) 

where c 2 = ci(E|Ci| r ) < oo; observe that the finiteness of the latter con- 
stant is due to (pTi"3]) . Finally, putting together (f2"7L5|) . ([277]), and Lemma [231 
in which one takes T = 1 and (3 = r/2 — 1, one obtains the theorem. □ 

From now on, our goal is to show that Relation (|2.5p holds, to this end 
one needs the following lemma. 

Lemma 2.5 Assume that to E [0,1] is arbitrary and fixed. Let {h n ) n is an 
arbitrary sequence of non vanishing real numbers which converges to and 
satisfies, for every n, to + h n G [0, 1]. Then, 

(i) for all fixed e > 0, conditional on the a -algebra Qx defined in 
the random variable 

__ Z(t + hn) ~ Z(tp) 
to ' hn |/j Jl/2+e 

has a centered Gaussian distribution with a variance given by, 

l-ma,x(t ,t +h n ) 

var[Rt , hn \g x ) = \h n \~ l ~ 2 ' / \$(X(s))\ 2 ds; 

J m'm(to,to+h n ) 
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(ii) moreover, under the assumptions that <I> does not vanish except on a 
Lebesgue negligible set and E\X(to)\ 2 > 0, one has, almost surely, 



l( ir{R to ,h n \Gx) > +00. 



Proof of Lemma 12. 5t Part (i) follows from (|2.9j) and the fact that condi- 
tionally on Qx, the stochastic process {^(i)}te[o,i] has a centered Gaussian 
distribution. Let us now prove that Part (ii) holds. The Mean Value The- 
orem, implies that, that there is a real number s n which belongs to the 
interval ( min(io, to + h n ), max(io, to + ^n)) an d satisfies, 



rma,x(t ,t +h n ) 

/ |$(X( S ))| 2 ds = \h n \\$(X(s n ))\ ; (2.16) 

Jmm(t ,t +h n ) 



1 1 2 

thus, using the almost sure continuity of the random function s 1— > \&(X(s))\ , 
one gets that, 

\hn\~ 1 / \$(X(s))\ 2 ds \$(X(t ))\ 2 . (2.17) 

Jmm{t,t+h n ) n->+oo 

Next, observe that the assumptions that $ vanishes only on Lebesgue neg- 
ligible set and that E|X(to)| 2 > 0, imply that, almost surely, 

\HX(t ))\ 2 >'0. (2.18) 
Finally it follows from (j2T7l) and (pTTB]) that 

\K\- l - 2t / |*(*(«))| ds +00. 



m'm(to,to+h n ) 



□ 



Now, we are in position to prove Relation (|2.5|) . 
Proof of Relation (12.51) : Our proof is inspired by that of Proposition 2.4 
in [3]. It consists in showing that for all fixed e > 0, there exists a sequence 
( r fc)fceN °f non vanishing real numbers converging to which satisfies: 

+ 1 r m 11 f 11, a \ Z (to + rk) - Z(t )\ a.s. fn-tn\ 
to + rk G [0, 1J tor all k, and : — -7— > +00. (2.19) 

\ r k ' fc^+oo 

To this end, it is sufficient to prove that there exists a sequence (h n ) n of non 
vanishing real numbers converging to which satisfies: 

t + K £ [0, 1] for all n, and — — = ^"J — — - — — ■> 0, 

R U h h n \Z(t + h n ) - Z(t )\ 

(2.20) 
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where the notation " — > 0" means that the convergence to holds in prob- 
ability. Indeed, assuming that (|2.20p is satisfied, then one can extract from 
(h n )n a subsequence denoted by (rk)kez + such that one has (|2,19p . Let us 
now prove (|2.20p . Denote by {h n ) n an arbitrary sequence of non vanishing 
real numbers converging to and such that to + h n € [0, 1] for all n. Observe 
that, one has, for all real number ij > 0, 

Pf <77|=lfl/ , \1 = E( E(l/ , \ Gx)) 



0.™n 

El 1(1 



Rt ,hn>Vv) 



Gx 



(2.21) 

Moreover, Lemma 12.51 entails that, almost surely, 



El/ \ 



+oo 



Gx) =\~ / exp(-x72)dx 
1 V vr.Vi(varK, h Jfo))- 1/2 



(2.22) 

and 

var (R t0 ,hJGx) +oo; (2.23) 



Thus, combining (|2.22p with (|2.23|) one obtains that, 



E l/ ,\ 

{Rt ,h n >Vv) 



Gx ) J- e" x2/2 dx = 1. (2.24) 







Finally, in view of (|2.2ip and (j2.24p . using the dominated convergence the- 
orem, it follows that (|2T2Dj) holds. □ 

3 Random series representation of {Z(t)} te ^ ^ via 
the Haar basis 

In order to state the main result of this section, we need to introduce some 
notations. 

• We denote by L 2 ([0, 1]) the usual Lebesgue Hilbert of the square inte- 
grable real- valued deterministic functions over [0,1]. 

• The Haar orthonormal basis of L 2 ([0, 1]) (see for example [19 1 112 } [23]). 
is the sequence of the functions: 

{<Po,o,-4>j,k : J £ Z + and k G {0, . . . , 2° - 1}}, 

defined as 

Vo.o = l[o,i) (3-1) 
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and 

ip jjk = 2 j / 2 (t, jL 2k+ i, - 1 [M1 w-iO, (3.2) 
where I5 is the indicator function of an arbitrary set S. 

• We denote by (fi,.F, P) the underlying probability space, that is the 
probability space on which the processes {-^(s)} s g[o,il > {^( s )}se[o,i] 
and {Z(t)} t £\ 0> i] are defined; moreover, for the sake of simplicity, we 
assume that this space has been chosen in such a way that: for all uj G 
n, Relation JT2J holds and G C"?([0,1]) for each 7 G [0,1/2). 
Also, we denote by 

{*o,o, hk ■ j G Z + and fc G {0, . . . , 2 j - 1}}, 

the sequence of standard independent Gaussian random variables de- 
fined, on this space as, 

So,o= f <po, (8)dW(8) = W(l) (3.3) 
J 

and 

A i>fe = /V J>fc ( S )d^( S ) = -2^(^(^)-2^(^±i) + Ty(A)). 

(3.4) 

Observe that, similarly to Lemma 2 in [3], one can show that there 
exist fi* C fl an event of probability 1, and a positive random variable 
C* of finite moment of any order, such that one has for all all 
j G Z+ and all k G {0, ... , 2> - 1}, 

|A i)fc (o;)| < C*(u)y/T+j; (3.5) 

we note in passing that the proof of this important inequality, mainly 
relies on Borel-Cantelli Lemma. 

• {K(t, •s)}(f jS )e[o,i] 2 ' 1S the stochastic field defined for all (t, s, uj) G [0, l] 2 x 
Vt as, 

K(t,8,u) :=*(X(8,u))l m ( 8 y, (3.6) 
thus, the process {-Z(i)}te[o,ii defined in (jl.ip can be expressed as, 

Z(t) = / if(t,s)dW(a). (3.7) 
Jo 

Observe that §3J^), dZHSD , dZH2D , imply that for all w G SI, 

sup |if(t,s,o;)| = sup \$(X(s,L}))\ < Ci(uj) < 00, (3.8) 
(t,s)e[o,i] 2 «e[o,l] 

where the last inequality results from the continuity of the functions 
s 1 y X(s,oj) and x h-> $(x). 
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• We denote by {&o,o(£)}te[o,ll the stochastic process defined for all (i, to) G 
[0, 1] x Q, as, 

b Qfi (t,u})= / K(t, s,Lo)tp ,o(s) ds; (3.9) 
Jo 

moreover, for all j G Z + and A; G {0, . . . , 2 J — 1}, we denote by 
{oj,fc(i)}t£[o,i] the stochastic process defined for all (t,uj) G [0, 1] x S7, 
as, 

a j:k (t,u) = / if(i, s,to)ipj i k(s)ds. (3.10) 

JO 

Observe that it follows from ([32]), (l3lH . (13^1 . (l3TTj) and ([32]), that 
for all to G 12 and for all t±, t 2 G [0, 1], 

|6 ,o(ii ! w)-6 ,o(i2,w)| <Ci(w)|ti-t 2 | (3.11) 

and for all j G Z+ and k G {0, . . . , 2 J ' - 1}, 

\a j>k (t u u) -a jik (t 2 ,oo)\ < C 1 {to)2i/ 2 \t l -t 2 \. (3.12) 



Now we are in position to state the main result of this section: 

Theorem 3.1 Let 7 G [0, 1/2) be arbitrary and fixed. For all J G Z + and 
(t,to) G [0,1] x $7, one sets, 

J 27-1 

Zj(t, w) = 60,0 (t, w)£ 0) oM + ^ ^ o;)A iife (a;). (3.13) 

3=0 k=0 

In view of 113.11]) and 113.12]) the trajectories of the process {^j(i)}*G[o,i] be- 
long to the Holder space C 7 ([0, 1]) since they are in fact Lipschitz functions; 
moreover, there exist fig an ^vent of probability 1 and a positive random 
variable D of finite moment of any order, such that one has for all to G Sl\ 
and J G Z + , 

\\Z(-,to) - Z J (.,^)|| C7([0il]) < 0( w )2-^(V2-7,«-l/2) ^TTJ, (3.14) 

where \\ ■ ||ct([o,i]) * s ^ e ^swa/ norm on C 7 ([0, 1]) (see Definition \2.1]) and 
where a G (1/2, 1] has been introduced in hl.2\) . 

In order to prove the latter theorem we need some preliminary results. The 
following lemma is a weak version of Theorem 13.11 

Lemma 3.2 We use the same notations as in Theorem \3.1\ Let t G [0,1] 
be arbitrary and fixed. When J goes to +00, the random variable Zj(t) 
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converges to the random variable Z(t) in the Hilbert space L 2 ($7), namely, 
one has, 

lim E(Z(t) - Zj{t)) 2 = 0. (3.15) 

J->+oo 

As a straightforward consequence, there exist O J t an event of probability 1 
included in SI* (recall that $7* is the event of probability 1 on which \3. 5\) 
holds) and a subsequence n \- > J n (a priori depending on t) such that for all 
wefijj, one has, 

lim Z Jn (t,uj) = Z(t,u). (3.16) 

Proof of Lemma 13. 2t Let K be as in (|3.6|) . Observe that in view of (|3.8|) . 
for all fixed (t,u) E [0, 1] x f2, the function K(t, :sh> K(t, s,u) belongs 
to L 2 ([0, 1]; ds). By expanding the latter function on the Haar basis, one 
obtains that, 

K(t,-,u) = bofl(t,u)tpo fi (-) + ^2^2 a i,fc(*5 w )V'i,A( , )> ( 3 - 17 ) 

j=0 k=0 

where the coefficients bo t o(t,u}) and aj^{t,uj) have been defined respectively 
in ()3.9p and (|3.10p . A priori the series in (|3.17p is convergent in the L 2 ([0, 1]; ds) 
norm, namely, 

lim / \K(t,s,u) - Kj(t,s,oo)\ 2 ds = 0, (3.18) 
where, for each J E Z+, Kj{t,-,ui) is the partial sum defined, for all s E [0, 1], 



Kj(t, s, w) = 6o,o(*) w)y2 ,o(s) + $Z a J.fc(*' w )V ; j,fc(s)- (3.19) 

j=0 fc=0 

Let us show that this series, is also convergent in the £ 2 ([0, 1] x 0; ds IP) 
norm, that is, 

lim e( f \K(t,s) -Kj{t,s)\ 2 ds] = 0. (3.20) 
J-H-oo Vio / 

In order to show that (j3.20p holds, we will use the dominated convergence 
theorem. It follows from (|3.17p . (|3.19p . Parseval formula and (|3.8p . that for 
all (oj, J) E Q x Z+, 

/•l 

|i^(t,s,w) -ifj(i,s,w)| 2 ds < / |^T(t,s,w)| 2 ds < C 2 (u;); 



thus, in view of (|3.18p and (|2.13j) . we are allowed to derive (|3.20p by making 
use of the dominated convergence theorem. 
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(J \K(t,s)-Kj(t,s)\ 2 ds) ; 



Finally, observe that, §31$, (pOjl . (|3T4"p . (|3~T3]l . (|3TT9j) . (fPjh and the 

isometry property of Wiener integral, entail that, one has almost surely, for 
all J G Z+, 

E ((Z(t) - Zj(t)) 2 |g x ) = J -i^j(t, S )| 2 d S , 
which in turn implies that, 

E(Z(t)-Zj(t)) 2 = E(E((Z(t)-Zj(t)y\g x ] 
= E 

therefore, using (|3.20p one gets (|3.15p . □ 

The following lemma provides sharp estimates for a,j t k(t). 

Lemma 3.3 There is a positive random variable A of finite moment of any 
order, such that for all U) E CI, t & [0, 1], j G Z + and k G {0, . . . , 2 3 — 1}, 

(i) when t > (k + 1)/2 J , one has 

\a j:k (t, u) | < A{uj)2-^ a+1 ^ , (3.21) 
where a has been introduced in il.ty) ; 

(ii) when k/2 j < t < (k + 1)/2 J , one has 

\a jtk (t,u})\ < A(uj)2 j/2 (t - k/2 j ); (3.22) 

(Hi) when t < k/2 3 , one has, 

a jyk (t,uj) = (3.23) 



Proof of Lemma 13. 3t First observe that (|3.23p easily results from 
(|3.6p and (|3.10p . Let us now prove that (|3,2ip holds, so we assume that 
t>(k + l)/2>. For all (y,u) G [0, 1] x O, we set 

X { ~ 1) (y,uj)= f\(X(s,oj))ds. (3.24) 



Then, using ([52]) . (|33|) . (j37[0]) and (pT24"|) . it follows that a jik (t,u) can be 
expressed as an increment of order 2 of the function 2^ 2 X^ ^ (•, uj): namely, 
when t > (k + 1)/2 J , one has, 



a jyk {t,uj) (3.25) 

-^(xr'(^±i, W )-At"(^.-)). 
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Next, applying the Mean Value Theorem to the function 

on the interval [0, 2~i~ ], it follows from (|3.25|) and (|3.24|) . that there exists 
z G (O^-^ 1 ) such that 

a j>k (t,u) = (*(X(^±± - z,J)) - *(*(^ - z, a,))) ; 

(3.26) 

then, applying the same theorem to the function x t— >■ $(x), on the compact 
interval whose endpoints are X ( ^j+l —z,ui^ and x(^~- — z, u^J , one obtains, 
in view of (pT26|) and (HTTP]) , that, 



X — — Z, U)\ — A — — : Z.LJ 



2i +1 'J V 2i 



(3.27) 



where 

C 2 (w) :=sup{|$'(s)| : \x\ < ||^||oo(^)}. 

Observe that the continuity of the functions s t— >■ X(s,u) and x t— >■ &'(x), 
implies that C% is a finite random variable, moreover (jl.3p and (|2.1ip entail 
that all its moments are finite as well. Next, combining ()3.27p with (|1.2p . 
one gets (|3.2ip . Let us now show that (|3.22p holds. It follows from (|3.2p . 
(J3SD and (I3TH that 



a 



it k(t,u>)\ < 2^ 2 / ds < d(u)^ 2 (t - fc/2 J ), 



where the random variable C\ has been defined in (j2.12j) : thus we obtain 

Lemma 3.4 Let 5 &e an arbitrary fixed positive real number. There is a 
constant c > 0, only depending on 5, such that for all x,y G (0,1], satisfying 
x < y, one has 



x 5 ^l + log^x- 1 ) < c/ Vl + log 2 (^ 1 )- 

Proof of Lemma 13. 4i The derivative over (0, 1] of the positive function 
5 i y s^ y/l + log 2 (s _1 ) is equal to 



log 2/ V lo g 2 21 °g 2 

Let us set ss = 2e~^ 1 . When ss > 1, then the function s H> \Jl + log 2 (s -1 } 
is increasing on (0,1]; thus taking c = 1, it follows that the lemma holds. 
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When s$ £ (0,1), then the latter function is increasing on (0, sg] and de- 
creasing on \ss 1] ; thus setting 



_ max se(0|1] {s 5 y 7 ! + log 2 (.s r )} 
nun a6 p 4) i] {s 5 V 1 + 1 og 2 (s -1 )} ' 

one obtains the lemma. □ 

Now we are in position to prove Theorem 13.11 
Proof of Theorem I3.lt First notice that it is sufficient to show that for 
all to £ f2* (recall that f2* is the event of probability 1 on which (|3,5p holds), 
J G Z + and Q S Z + \ {0} one has, 

\\Z J+Q (.,«,) - Zj(;t0)\\ CJm]) < £,( fa ,) 2 -Jniin(l/2- 7 ,a-l/2) (3 . 28) 

Indeed, (|3.28[) implies that (Zj(-, w)) j gZ+ is a Cauchy sequence in the Ba- 
nach space (7 7 ([0, 1]) and, as consequence, it converges, in this space, to 
some limit denoted by Z(-,to). Next let Q 2 be the event of probability 1 
defined as 

^= n n iv 

ge[o,i]n<Q 

where fij C fi* is the event fi* t introduced in Lemma [3.21 when t = q. Thus 
it follows from the latter lemma, that for each to £ ^ an d au Q ^ [0, 1] fl Q 
one has Z(q,co) = Z(q,to); this is equivalent to Z(-,to) = Z(-,to), since 
Z(-,to) and are continuous function. Then, letting in (|3.28[> . Q goes 

to +oo while J fixed, one gets (|3.14p . 

From now on, our goal will be to show that (|3,28p holds. To this end, in 
view of (|2.ip . it is sufficient to prove that there exist two positive random 
variables D\ and D 2 of finite moment of any order, such that one has for all 
wG!l*,JeZ + and Q £ Z + \ {0}, 

\Zj+Q(h,u) - Zj(ti,uj) - Z J+Q (t 2 ,<o) + Zj(t 2 ,(o)\ 

sup J : : 

o<ti<t 2 <l \t\-t 2 p (3.29) 

< jDl ( w ) 2 - Jmin ( 1 /2-7,«-l/2) ^Jl^Tj 

and 

sup \Z J+Q (t,to)-Zj(t,to)\ < ZM^)2- Jmin(1/2 - 7 ' Q - 1/2) ^TTJ. (3.30) 

0<t<l 

We will only show that (|3.29p is satisfied, since (|3.30p can be obtained in the 
same way. Let t\ < t 2 be two arbitrary and fixed real numbers belonging 
to the interval [0,1]. Observe that in view of (|3.13|) . using the triangle 
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inequality, one has, 



\Zj + Q(ti,u) - Zj(ti,u) - Z J+Q (t 2 ,u) + Zj(t 2 ,uj) 
J+Q 2^-1 

= | ^ K'.fctA'^) ~ a i,fc(i2,w))A :) - fc (w) 

j=J+l k=0 

+00 2^-1 

j=J+l k=0 



(3.31) 



Let us now give appropriate bounds for 

2J-1 

Uj(t\,tv,u}) = j ^ (a j>k (ti,oj) - a jtk (t2,u))^j,k(u) (3.32) 

fc=0 

and derive form them (|3.29p . We denote by jo G Z + the unique integer such 
that 

2- j °- 1 < \h - t 2 \ < 2~ jo . (3.33) 

Also for all t £ [0, 1] and j £ Z+, we denote by k(j,t) the unique integer 
belonging to {0, . . . , 2 J — 1} such that 



k(j,t) < t) + 1 



2J 



2-' 



with the convention that 



Observe that when t £ [0, 1), 



1) = V - 1. 



(3.34) 
(3.35) 

(3.36) 



[•] being the integer part function. Also, observe that ()3.32p . (|3.23p . and the 
fact that ajk(t,(J) does not t when t > (k + l)/2- 7 (see (|3.25p ). imply that 



fc(i,ts) 
k=k(j,h) 



(3.37) 



Let us now study the following two cases: j < jo and j > jo. First assume 
that j < jo, then (|3.33p implies that 2~ 3 > \ti — t 2 \ and, as a consequence 
that k(j,t 2 ) £ {k{jM)MjM) + l}. When, k(j,t 2 ) = fcfrti), it follows 
from flEgp , (|53D and (J3H2D, that 



Uj(ti,t 2 ,uj) 



,S(j,ti)( tl ' w ) " °i.*(j,t 1 )( t 2,w)||A i)£(j . tl) (a;)| 



< Gi(w)|ii- i 2 |2 J '/ 2 x/T+7, 
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where G\{oS) := C*{bS)Ci(oS). When, k(j,t2) = k{j,t\) + 1, putting together 
([37371) . (133]) . (j3~12l) . (13321) and the fact that (k(j,h) + l)/2-?' G [ti,t 2 ], one 
obtains that, 

Uj(ti,t 2 ,u)) 

^ K'kUM)^'^ ~ a jMi^^)\\hmt l) ^)\ + l^fco-, fl )+i( t 2^)||\fe( ij t 1 )+i( CJ )l 

< Gx{uS)\ti - t 2 \2^ 2 y/T+j + G 2 (w)(t2 - (fc(j, ti) + l)/2 J )|ti - t 2 |2^ 2 yiTj 

<G 3 Miti-t 2 |2^ 2 yrT7, 

where G 2 (w) := A(w)C*(a;) and £3(0;) := Gi(cj) + G 2 (co). Thus, we have 
shown that there is a positive random variable G4 of finite moment of any 
order, non depending on 7, jo, t\ and t 2 , such that for all j < jo, one has, 



Uj(t\,t 2 , 



\h-t 2 \^ 

Therefore, in view of (|3.33|) . for each integer J satisfying < J < jo, one 
has, 

jo TT , \ jo 

E 6 h :^ <G^)\ tl - t2 \^ e 2 ^0T7 

< ^(72 - i)^ 1 G 4 (^)|ti - tsi 1 - 7 2^/ 2 yrT7^ 

< V2(V2 - l)" 1 G 4 H|ti - t 2 | 1/2 - 7 Vl + log 2 (l*i - tal- 1 )- 

(3.38) 

Next, it follows from lEOB]) . the inequality 2~ J > 2~i° > \h-t 2 \ (see ([3331) ) 
and Lemma 13.41 (in which one takes 5 = 1/2 — 7, x = \ti — t 2 \ and y = 2~ J ), 
that 

E y i>t /^ ) <G5M2- J ^-7)VT+J, (3.39) 

where 

G 5 (w) := cV2(V2 - l) _1 G 4 (w), 

c being the constant introduced in Lemma 13.41 Let us now study the case 
where j > jo- Observe that in this case, in view of Relations (|3.33|) and 
(|3,34p . one necessarily has k{j,t\) < k(j,t 2 ) and 

2 j \h-t 2 \>l. (3.40) 

Also observe that, using (|3.35p . (|3.36j) and (|3.40p . one obtain that 

KiM) - Hj,h) < 2H 2 - 2H X + 1 < 2-? +1 |t 1 - t 2 \. (3.41) 
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It follows from (|332l) . (|3T25l) . Lemma ET3l (pOT]) . m and (IXT2D . 

that 

C/,-(ti,*2,w) < |a J -,fcy )tl )(*i>w)||A J . i j E(j . jtl) (w)|+ 2 | (*2, 1 1 Aj ifc (w) | 

k=k(j,ti) 

< 2G 2 (co) (2^/2 yfj- + | tj _ t2 \ 2 -i(a-i/2) 

(3.42) 

Thus (|3.42p and (|3.40p imply that for all j > jo, one has, 

< 2G2M / 2 _ i(1/2 _ 7) ^^ +|ti _ i2|1 _ 72 _ i(a _ 1/2) ^ T ^ 



|*i -t 2 | 7 

< G 6 (w)2--'' min ( 1 / 2 -^ Q - 1 / 2 ) v/l+7- (3-43) 

where := 4G2(w). Next, observe that one has for J £ Z + , 

+00 

^ 2 -i min ( 1 /2-7,a-l/2) yiTj 



2 -Jmin(l/2- 7 ,Q-l/2) ^[T^J ^ 2 

i=J+l 



00 

(j-J) min(l/2-7,a-l/2) 



i + i 



1 + J 

< Cl2 - Jmin ( 1 /2-7,«-i/2) 7T^7j ; (3.44) 
where the constant 

+00 

Cl := J- 2 -imin(l/2- 7 ,a-l/2) < +QO _ 

1=1 

Thus combining ()3.43p with ()3.44j) . it follows that for all J > jo, one has, 

g Uj(tut 2 ,u) ^ G7( ^ 2 -J m m(l/2- 7 ,a- 1 /2) ^/j^ (3 45) 

where Gj(u) := ciGg(w). Let us now show that, for all J G Z+, one has, 

V Mllhl^ < g (cj)2 -Jmin(l/2- 7 ,a-l/2) ,3^ 

where -Di(cj) := (1 + 0)^7(0;) + G$(ui), c is the constant introduced in 
Lemma l3~4l and G^{ui) has been introduced in ()3.39p . It is clear that (|3.45p 
implies that (|3.46p holds when J > jo, so from now on, we assume that 
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jo > 1 an d that J is an arbitrary nonnegative integer satisfying J < jo- It 
follows from ()3.45p that 

g Mhhl^l < G?(fa;)2 - 3 - min(l/2- 7 ,a-l/2) 

Then using Lemma 13.41 (in which one takes 5 = min(l/2 — 7, a — 1/2), 
x = 2~i° and y = 2~ J ), one obtains that 

Y ^fr'* 2 '") < cG (a;)2 -/min(l/2- 7 ,a-l/2) yTTjTj^ (3.47) 
.41. t 1 -t 2 T ~ 

Next combining (|3.47|) with (|3.39p . it follows that (|3.46p holds in the case 
where J < j . Finally (pUiil . (pU?2l and (EnBj) imply that (ET29j) is satisfied. 
□ 



4 Simulation of {Z(t)} ie [ ,i] v * a the Haar multires- 
olution analysis 

Our algorithm for simulating {Z(t)} te [$ u mainly relies on Theorem 13.11 
which allows to approximate {^(t)}ie[o,i] ) f° r <^ large enough, by the pro- 
cess {Zj(t)} te \ 0t i\ defined in f)3. 13|) as a finite of sum. First we will give an 
expression for the latter process which makes it rather easy to simulate, to 
this end, we need to introduce some notations. 

• For all fixed integers J £ Z + and I G {0, . . . , 2 J — 1}, the function ipjj 
is defined as, 

ifj ) i = 2 J lH,^ l±ly (4.1) 

• For all fixed integer J G Z+, we denote by {Sjj : I G {0, . . . , 2 — 1}} 
the finite sequence of independent standard Gaussian random variables 
defined as, 

hi ■= f Q <PjM dW(s) = 2 J / 2 (W( l -±1) - W(±j)) . (4.2) 

• For all fixed integers J G Z + and /£{0,...,2 J — 1}, the stochastic 
process {frj,z(^)}te[o,i] IS defined for all (t,u) G [0,1] x Q, as, 

bj,i(t,uj)= K(t,s,w)v9 Jj/ (s)ds, (4.3) 
J o 

where K(t,s,uj) has been introduced in (|3.6|) . 
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The following proposition provides a nice expression for the process {Zj(t)} tE \ 0t x\ . 

Proposition 4.1 For all fixed integer J 6 Z +; let {Zj(t)} t€ [ ^ be the 
stochastic process defined in i3.13\) . Then one has for all t E [0,1], almost 
surely, 

Zj{t)= b w {t)5 w . (4.4) 

1=0 

Observe that, Relation (j4.4[) . also holds almost surely for all J G Z + and 
t £ [0, 1] (i.e. this relation is satisfied on an event of probability 1 which does 
not depend on J and i), since the trajectories of the processes {Zj(t)} te [ ^ 

and {Xw=o _1 bj+ij(t)5j+i t i } , are with probability 1, continuous func- 

tions. 

Proof of Proposition I4.lt The main ingredient of the proof is the Haar 
multiresolution analysis (see for example |19[ [T2"| 125] ) of the Hilbert space 
L 2 ([0, 1]); that is the increasing (in the sense of the inclusion) sequence 
(Vj)jzz + of the finite dimensional subspaces of L 2 ([0, 1]) defined as, 

V := Span {990,0} 

and for all integer J > 1, 

Vj := Span {ipo fi , tpj,k ■ 3 € {0, . . . , J - 1} and k £ {0, . . . , 2 j - 1}}, 

where the orthonormal functions c^o,o an d ipj t k have been introduced respec- 
tively in (|3.ip and in (j3.2j) . Relations (|3.17p and ()3.19p imply that for every 
fixed J £ Z + and (t,oj) £ [0, 1] x the function Kj{t, -,w), can be viewed 
as the orthogonal projection (in the sense of the usual inner product of 
L 2 ([0, 1])) of the function K(t, -,uj), on the space Vj + \. On the other hand, 
it is known (see for example [T9J [121 E3]) that, the finite sequence, 

{<p J+1>l , Ze{o,...,2 J+1 -i}}, 

forms an orthonormal basis of Vj+i. Therefore, in view of (|4.3p . one has, 
for all s e [0,1], 

2 J + L -1 

Kj(t,s,u)= b J+h i(t,uj)ipj +h i(s). (4.5) 

1=0 

Finally, using the fact that, 

Zj(t) = f Kj(t,s)dW(s) 
Jo 
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as well as Relations (ETKJ]) . (ET41 . (|43|) and (g2D, one obtains (|Q|> . □ 



Now, we are ready to describe the main step of our algorithm for simu- 
lating {Z(t)}i 6[0il] . 

Main steps of our algorithm for simulating {Z (t)}te\o,i\ : 

(1) We take J large enough and we simulate the finite sequence 

{Sjy. /G{0,...,2 J -1}} 
of the standard independent Gaussian random variables defined in (|4.2p . 

(2) We simulate 

•>2J 



observe that in the case where the centered Gaussian process {-X"(t)}te[o,i] 
is a multifractional Brownian motion; such a simulation can be done 
by making use of one of the efficient methods described in [5] . 

(3) Noticing that, for all I G {0, . . . , 2 J — 1}, and m G {l + l,..., 2 J }, 
Relations ([43]) . (|43j) and (J3ISD imply that, 

i+i 



b 4^)= 2J/2 I7 ^ x ^^ 

2J 



we approximate the latter integral, by the Riemann sum 

9=0 V 7 



(4.6) 



On the other hand, observe that (|4.1|K (|4.3p and (|3.6() entail that for 
each m G {0, ...,/}, 

m 1 



(4) Thus, in view of (|4.4p . for all m G {1, . . . , 2^}, we approximate -Zj-i (^7) 
by 



Mf) -EMfK < 4 - 7 ' 

Z=0 



Then we simulate 
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by using the fact that 

Zj-iQj) =h,o(±j)6j,o (4.8) 

and the induction relation, for all m G {2, . . . , 2 J }, 

~ f m\ -~- /m—l\ -> ( m\ „ 

= Zj- l (—)+b J ^ 1 (- ¥ )5 J , m - 1 . (4.9) 

Observe that (gSJl and (02]) easily result from (02]) and flUj) . 
(5) Finally, by interpolating the 2 J + 1 points 

(0,0); Qp^-i(^)) ; ... ; ; (1.^-1(1)) 

we obtain a stochastic process {Zj_i(t)} ie j Q ^ which, in view of The- 
orem [37T1 satisfies the following property: for all fixed 7 G [0,1/2), 
there exists a random variable D' of finite moment of any order, such 
that one has, almost surely, for all J big enough, 

\\Z - £j-i|| C7([0)1]) < £) / 2- J,nin ( 1 /2-7.«-i/2) VT+J, 
where || • [|c7([q,i]) is the usual norm on the Holder space C 7 ([0, 1]). 

□ 
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The Hrust functional parameter H(t)=0.6+0.2t 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

time 

The mBm with Hurst functional parameter H(t)=0.6+0.2t 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

time 



Figure 1: Graph of the function H(-) such that H(s) = 0.6 + 0.2s for all 
s G [0,1], and a simulation of a trajectory of a multifr actional Brownian 
motion {^(s)} s e[o,i] with Hurst functional parameter H(-). 
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The Irajctory of the MSVM, Phi(s)=0.5+0.5s 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

lime 

The trajectory of the MSVM, Phi(s)=Sin{s) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

time 



Figure 2: Simulation, in the case where <J?(x) = 0.5 + 0.5a; for all x £ R, 
and also in the case where = sin(x) for all x G R, of a trajectory of 

the log price process {^(i)}te[o,l] generated through (jl.ip by the multifrac- 
tional Brownian motion {-X"(s)} s e[o,i] m Figure [IJ The shapes of these two 
simulated trajectories tend to confirm the fact that the pointwise Holder 
exponent of {^(i)}te(o,l] does not change from one place to another and is 
equal to 1/2 (see Theorem 12. 2p . 
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